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Abstract. This paper investigates the evolution of collapsing FRW models with a scalar 
field having the potential which arises in the conformal frame of high order gravity theories, 
coupled to matter described by a perfect fluid with energy density p and pressure p, obeying 
a barotropic equation of state. The solutions are shown to evolve generically to a singular 
state in a finite time and they are used as sources for radiating objects undergoing complete 
[ — ■ gravitational collapse. Although these singularities may be naked in some special case, it is 

\ shown that generically a black hole forms. 
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1. Introduction 

The study of Friedmann-Robertson-Walker (FRW) solutions has always been a central 
topic in relativistic cosmology. In this sense, a great deal of attention is paid to investiga- 
tion of global structures and, waiting for a satisfying theory of quantum gravity, spacetime 
\ boundaries may be described in classical terms only, studying existence of singularities when 
gravity is coupled to different kind of matter tensors. Especially dealing with early universe, 
a huge literature exists on scalar field cosmologies, both in the free case [01], and in the more 



complicated situation when a potential is added (see for instance [g, [12], [jj]) and initial (i.e. 
big-bang) singularity existence is investigated. 

Scalar field solutions are of great importance also in relativistic astrophysics - massless 
solution is a useful model in gravitational collapse studies because the evolution equation 
in absence of gravity is free of singular solution. On the other side, adding a potential to 
the model may result again in singularity formation, which here represents the endstate of 
collapse JT3|]. In this context, the gravitational collapse of self-interacting homogeneous 
scalar field models has been recently analyzed in [TTT]], where it is proven that for a wide class 
of potential the evolution is generically (i.e. up to a zero-measured initial data set) divergent 
in a finite comoving time. 

In this paper flat FRW spacetime is considered, where scalar field is nonminimally cou- 
pled to ordinary matter described by a barotropic fluid. Motivations come from higher order 
gravity theories (HOG) derived from Lagrangians of the form 

(1.1) L = f(R) 
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where / is an arbitrary smooth function and L m is the matter Lagrangian depending on 
the matter fields Nonlinear theories of gravitation are a well-established field of research 



since the pioneering works [[H IT5L 21 1; see also P] for a thorough analysis of interpretation 



issues. It is well known that under the conformal transformation 

(1.2) g»v = f{R)9nu, 

the field equations reduce to the Einstein field equations with a scalar field as an additional 
matter source, namely 

(1.3) G IH , = T ia ,{g,<j>) + T IH ,{g,V), 
where 

T, u (g, 0) = d^dA - l -g„ v [{dct>f - IV (0)] , 

and 

(1-4) 0=y|ln/'(i2), 
allowing to write the potential of the scalar field as 

V(R(cf>)) = j^— 2 (Rf'-f). 

From now on let us work exclusively in the conformal frame and for simplicity drop the 
tilde from all quantities. Ordinary matter is described by a perfect fluid with equation of state 
p — (7 — l)p (with < 7 < 2) and field equations reduce to the following system: 



(1.5b) H = - -<p 2 -^rp + 



1 -,2 7 k 
2* ~2 P+ ^ 



(1.5c) + 3H<j> + V'{<f>) = ^yl 2 ^ 

4-37 • 

(1.5d) p + 3 1P H = -j=±p<P 

(see for example |16|]). Here H = a/ a, where a (i) is the scale factor of FRW and overdot 
denotes differentiation with respect to time t. It is worth stressing here once for all that this 
approach is different from the one used in previous investigations [0, [17p where the vacuum 
Lagrangian is conformally transformed into the Einstein frame and the matter Lagrangian is 
added (see [|8|] for a comparative discussion about the different approaches). 

Of course, for both cosmological and astrophysical purposes, it is important to study the 
solutions of fll.5aj )-( |l3dl ) relaxing as much as possible hypotheses on the potential, although 



exponential potentials, for instance, have been extensively treated in literature. In the present 
paper, the class of potentials studied is described in technical terms by Definition [23| . Poly- 
nomial potentials with even leading terms, exponential potentials - with small logarithmic 
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derivative value - and potentials arising in conformally related theories of gravity, such as 

HH 

(1.6) V ((/>) = Vo(l-e-V^^ 2 



belong to the class considered. In Section 0, casting the problem into a dynamical systems' 
framework, it is proven that the qualitative behavior of the solution is similar to the case 
without matter, i.e. solutions completely collapse in finite time up to a zero-measured initial 
data set (Theorem Responsible for the formation of this singularity is the scalar field, 
which dominates over the energy density of the fluid even when the latter diverges. 

In Section |3| the solutions are used as source to build collapsing object models where the 
exterior is given by a generalized Vaidya solution. The endstate of these models is studied for 
such objects, finding that the collapse generically ends into a black hole. We stress the fact 
that naked singularities may arise from these models, choosing initial data in the non generical 
set excluded by Theorem |2T2|. In this sense, the result is qualitatively similar to what shown 
in [ |TT] ] - for the class of potentials considered, weak cosmic censorship hypothesis is not 
affected when a barotropic perfect fluid is added to the model. 

2. Collapsing solutions in the flat case 
In the following, we are going to study the future late time behavior of solutions of ( |1.5a| )- 



( |1.5dp , with k — 0, such that they are collapsing at initial time of observation, i.e. H(0) < 
0. Since 7 is non negative, ( |1.5bD implies H(t) < 0, Vi > 0. Regardless of the physical 
derivation of the potential V ((/)), discussed in the previous section, the results we are going to 
state will hold for the following general class of potential functions V(0). 

Definition 2.1. We say that V belongs to the set € if the following conditions are satisfied: 

(Al) (Structure of the set of critical points) The critical points of V are isolated. 

(A2) (Existence of a suitable bounded "sub-level" set) There exists a < b e R such that 

> b V'(<f>) > 0, < a V'(<f>) < 0. 

Moreover, setting V* = max^j] V, it holds 

lim V{<j>) > K, lim V{<t>) > K- 

</>—>— 00 </>— »+oo 

(A3) (Growth condition) The function defined in terms of V and its first derivative V as 
satisfies 

(2.2) 3 lim \u{(f>)\ < 1, 

(j>— >±oo 

(2.3) 3 lim u'(0)(=O). 

The proof will be carried out translating the original problem into a compact framework, 
in order to apply classical results of dynamical systems. This approach, already used in 
[Q] - where gravity was minimally coupled to a scalar field under the action of a potential 
only - seems more promising in this framework, where a perfect fluid is added, than the 
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one used in [|TT|], where everything was reduced to the study of a single second order ODE. 
The present approach, however, makes necessary a further assumption on V(<f>) allowing the 
compactification of the problem, that will be performed once Theorem pj] below is proved. 



Assumption 2.2. We assume that for some constant M > b, there exist a C 2 map f((f)) 
[M, +00 [^]0, s ] such that f is monotone decreasing (and then \im ( p^ +00 f(<p) = 0, f(M) - 
so), there exists lim^ +00 /'(</>) (= 0) and moreover 

u'U) 

(2 - 4) , lim 77tS = °> 

0_H.BC f'((f)) 

(2.5) 3 lim ^}g1. 

0-+oo /' (0) 

Definition 2.3. We say that a potential V E £ belongs to the class £ if assumption 
holds for V((f)) together with a similar assumption for the behavior at negative infinity, i.e. 
supposing also the existence of a function (that we will call / again) / :] — 00, N] — > [— s , 0[, 
where N < a, with similar properties as in Assumption [2.2| opportunely adapted to this case. 



Remark 2.4. Although the class £ contains many relevant examples in cosmology such as 
polynomials with even leading term, exponential functions e A * (with A < y/E, to satisfy 
growth condition (A3)), or even the potential fll.6D arising in the conformal frame theory 
discussed in the introduction, one can exhibit examples showing that <£ \ £ is non empty: 
take for instance V(4>) = e <t> Sl (<t>)+ cos $ as <\> > b, which]] satisfies the growth assumption (A3) 
but does not admit a function /(</>) as in Assumption ^E2\ 

Although the main Theorem ^2] stated in Subsection ^2] will require V E £o> for some of 
the preliminary results (for instance, Theorem |2TTD it will suffice V E C 

Let us start the study of the problem introducing the unknown functions 
(2.6) * = -, y=-, z = — , 



and the time variable change dr = —Hdt. The system ( |1.5a| )H ]T3d] ) for k = becomes 

(2.7a) * = -y, 

dr 

(2.7b) ^ = _I a; ( 2/ 2 + 7Z 2 )5 

(2.7c) ^ = V'^)x 2 + 3y - \y* - z 2 {a + ly), 

(2.7d) ^ = ~z [y 2 -ay + 7(2^ - 3 

where we denote a := -7^-- Equation ( |1 .5a| ) gives the constraint 

(2.8) V{(t>)x 2 + l -y 2 + z 2 = ?>, 

which is invariant by the flow of the above system. We also observe, from ( |2.7bD , that the set 



{x = 0} is invariant by the flow and so sign(x) is also invariant. This fact guarantees that for 



'The function Si(4>) is the sine integral function, defined as J Q sind/fl d$. 
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collapsing solutions corresponding to curves with x < 0, x(r) will remain non positive for 
all r > 0. Let I C [0, +oo) be the maximal right interval of definition. The following crucial 
result holds. 



Theorem 2.1. Let V E <£. For almost every solution 7 = (<f), x, y, z) of ( |2.7aD -( p7d| ), swc/i 



z7ia? (2.8) and x < are satisfied at r = (anJ zTzen Vr E I J, zY /lo/ds 



limsup |0(t)| = +00. 

t— >sup I 

Proof. By contradiction, let 0(t) be bounded. Since x is negative and monotone, it is bounded 
on I and so, using (|2~%|) and the fact that V(<f>) is bounded from below (possibly by a negative 



value), the curve 7 is bounded, which means that sup I = +00. Let a>(7) be the omega-limit 
set of 7 which Jl8| , Section 3.2] is a nonempty, connected and compact set invariant by the 
flow of fl2.7a| )-( |277cl ). By continuity also ( |2.8| ) is satisfied on ^(7). 



Let p e u(7) and 7 P be the solution through p. Again, since x is monotone increasing on 
7, then x evaluated on the omega-limit set ^(7) must be constant, equal to the limit value of 
x on 7. This implies that x is constant on 7 p , which gives 

(1) either x = 

(2) or y = z = 0. 

Let us consider separately the two situations. In the first case, using (|2.8|) in (|2.7c[) we get 



(2.9) ^ = ^ 3 -^^l-^»- a 

and so y(r) (which is bounded) tends to an equilibrium point of (|2~13|). Equilibria are given 
by y = ±V6 and y = therefore, the only admissible case happens when = 0, 
that is 7 = 4/3, since otherwise y would tend to a nonzero value and, using fl2.7a| ), <p would 
be unbounded. In the admissible case 7 = 4/3 it is easy to see that 7 P reduces to the point 
x = y = 0, z = —V3 (we exclude the positive value by continuity). To prove that solutions 



of the system do not approach this equilibrium point generically, we use ( |2.8[ ) in ( |2.7bD - 



Q2.7cD and study the unconstrained system given by ( P.7aD -( j277cj ). It can be easily seen that 



eigenvalues of the equilibrium point ((f) = </> , x = 0, y = 0) are given by 0, —2 and 1, so 
there exists an unstable manifold at the equilibrium with strictly positive dimension. 

Then, it must be still analyzed case (||) where y = z = 0, which means <p = 0o constant 
and x = —-\/3/V(4>o). This fact implies V'(4>q) = and it is better to use the constraint 
( |1.5a| ) to study the behavior of solutions of 



"i(p + \v 2 + V(<P) ) r -V'(,,) + n f >. 



P = ~P I aw + 7^/3 (p+^v 2 + V(4>) 

in the unknown functions ((f), v, p) of t, near the equilibrium point (0 Q , 0, 0) (note that the 
system is regular at this point, since V((f) ) is strictly positive). The eigenvalues are given 
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by 7a/3K(0 o ) and \ ^3V((f) ) ± y/3V((j) ) - 41/"(0 O )J , which means that there exists at 

least one eigenvalues with positive real part and, again, there exists an unstable manifold at 
the equilibrium with strictly positive dimension. This completes the proof. □ 



Theorem |27T] suggests that one needs to find a suitable coordinate transformation, which 
maps a neighborhood of <p = oo onto a bounded set. Assumption |2.2| precisely performs this 
task, therefore from now on we will assume that the potential V(4>) belongs to the class £ 



(see Definition |23| ). 

Using the above hypothesis, we perform the change s = f(<p) and observe that = 
— = —yf'((fi(s)), where is the inverse function of s(4>) introduced in Assump- 

tion |2]2|. Observing that the above relation, together with ( |2.7d| ), does not contain x, we 



can eliminate x from ( |2.7cj ) using the constraint ( [2.8D and then we end up at the following 
unconstrained system: 



(2.10a) 

(2.10b) 
(2.10c) 



dy 
dr 

dz 1 r 
z\y 

ds 
d7 



-y 2 ) (y + V6u((P(s))) - z 2 {ly + a + V6u(4>(s)) 



ay + ^(z 2 — 3 

yf'(d>(s)). 



Because of assumption pT2| , we can extend the functions u((j)(s)) and /(0(s)) up to s = 0. 
Then, due to constraint ( [2.8D and to the fact that V(4>) is positive as — > oo, we will consider 
as admissible set for the system ( |2.10a| )-( [2.10c[ ), the set 



{(y,z,s) : \y 2 + z 2 <3, se[0,s ]}. 



The assumptions made on V((J)) allow us to consider the same problem ( [2. 10a| )— ( |2. 10b| ) 
on the set f2_ = {{y,z, s) : \y 2 + z 2 < 3, s 6 [— s ,0]}. Let us observe that the set 
{s = 0} may be seen as a subset of both fi + and but its meaning is different between the 
two situations, since it represents the flow at positive and negative infinity, respectively. To 
distinguish between the two cases, we will refer to them as the sets S + and XL, respectively. 

In this way we have reduced the study of our dynamical system in a compact set. To study 
the global properties of the flow, however, we need to take into account the behavior of the 
solution when the scalar field </> lies in the set [N, M). To this aim we first observe that, 



choosing initial data of ( |2.7aD -( ]2~77d| ) such that x(0) = x < 0, then x{r) 2 < x\, Vr > 0. 



Therefore we need to perform a C 2 junction of the sets Q + and Q_ described above, with a 
suitable compact subset f2 of {(y, z, 4>) : \y 2 + z 2 < 3 + Cx\, <fi G [N, M}}, where the 
constant — C (that can be chosen to be nonnegative) bounds V((f) from below, in such a way 
that the flow for t > completely lies in the set f2 = U fi U see Figure [[]. 



Definition 2.5. We say that a curve 7 : I — > f2 is a solution of the problem V if it solves 
( |2.10a| )-( |2.10c[ ) when 7 lives in f2 ± and it solves the system given by equation ( |2.7a| ), ( |2.7cp 
and Q2.7d| ) when 7 is in Vl , where x is given by the constraint (£ 



Remark 2.6. It is a simple consequence of the above construction, that every solution of the 



initial problem ( |1.5a| )-( [l3d] ) (with k = 0) such that H (0) < x , corresponds uniquely to a 
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Figure 1 . The set O is defined as the union of Q± with a suitable subset of 

{(y, z, 4>) : \y 2 + z 2 < 3 + Cx 2 Ql (p G [N, M)}. Note that coordinate chart on 
Q± and fi are given by (y, z, s) and (y, 2, 0), respectively. 



curve 7 in f2 solution of P. Of course, the converse may be false, because the set of solutions 
of V may also include curves lying in the set 

m = £+ U {(y, z, s)6fl + : \y 2 + z 2 = 3} U {(y, z, s) G f2_ : |y 2 + z 2 = 3}UL. 

The following lemma states a useful property of cu-limit points of the flow. 

Lemma 2.7. Let ^ &VLbe a generic solution of problem V and letp G u){"f) an uj-limit point 
for 7. Ifp G fi+ U fi_, p G Sfi. 



Proof. It must be shown that ^(7) H (0 + U C First, we observe that Theorem ^T] 
states the existence of a point q G ^(7) such that q G S. That means that lim t _^ sup i = 0, 
because the constraint ( |2.8[ ) is satisfied on 7 and then on ^(7). Therefore, if p G ^(7) H (0 + U 
f2_) is such that s ^ 0, then V is finite and again from ( |2.8| ) we get |?/ 2 + ^ 2 = 3, which 



implies p G □ 

The above lemma states, in rough words, that a generic solution of problem V approaches 
the set dVl. By the word "generic" we mean, in this context, a solution corresponding to a 
generic - up to zero-measured set - choice of initial data. There can be solutions such that the 
scalar (j) approaches, for instance, a local extremum of the potential V {(f)), but these solutions 
correspond to special choices of the initial data and are unstable with respect to this choices. 

2.1. The flow near 5f2. In order to study late time behavior of collapsing solutions, we need 
information on the structure of ^(7). Most of this information can be obtained studying the 
flow near 3f2. In this section we will list equilibria of the system fl2.10a| )-( |2.10cj ) and their 



properties. It will turn out that the phase portrait significantly changes, especially at S ± , 
depending on the values of 7 and w(0), the limit value of u(4>(s)) as s — > 0. The different 
situations depending on these parameters is summarized in Figure [ZJ each region corresponds 
to a particular phase portrait, sketched in Figure |3|. 

In the following, we will denote points of £± by (y, z, ± ) . 



8 



R. GIAMBO 




Figure 2. Summary of the different phase portrait of the flow near S±, ac- 
cording to the values of 7 and u^). The roman numbers correspond to situ- 
ations depicted in Figure [| 



(1) A±(-V6,0,0 ± ). The eigenvalues are Ai = -5+37, A 2 = 6(m(0 ± )-1), A 3 = 0, with 
corresponding eigenvectors vi(0, 1, 0), v 2 (l, 0, 0), v 3 (0, 0, 1). To find the behavior 
of the flow in the direction v 3 , note that ( |2.10c[ ) implies ^ = \/6 on v 3 , from which 
we get ps y/Er + O . This means that A + is attractive along v 3 and A- is repulsive 
along the same direction. Then we have the following three subcases, according to 
the value of 7: 

(a) 7 > 5/3. The point A± is unstable and the only stable direction tangent to £± is 
given by v 2 . 

(b) 7 = 5/3. The central manifold is 2-dimensional and it easy to see that 4 s = 
\/Qz 2 {\ — u(s)) > 0, so the point is repulsive along the direction Vi. 

(c) 7 < 5/3. The point is attractive along directions in the plane £± and so A + is a 
stable equilibrium point (whereas A- is unstable because repulsive along v 3 ). 

(2) i3±( v / 6,0,0 ± ). The eigenvalues are \ x = -1, A 2 = -6(1 + m(0 ± )), A 3 = 0, with 
corresponding eigenvectors Vi (0,1,0), v 2 (l,0,0), v 3 (0, 0, 1). Arguing as before, 
ps —y/6r + 0o along v 3 and so £L is a stable equilibrium point, whereas B + is 
unstable along the direction v 3 . 

(3) C±(- v / 6m(0 ± ),0,0 ± ). The eigenvalues are given by Ai = 3(1 - m(0 ± ) 2 ), A 2 = 
|(37 — 4^(0*) + 37m(0 ± ) — 6m(0 ± ) 2 ), A 3 = 0, with corresponding eigenvectors 
vx(0, 1, 0), v 2 (l, 0, 0), v 3 (0, 0, 1). The equilibrium point is unstable since Ai > 0, 
but it is useful to observe stability along directions on E± and so the sign of A 2 must 
be investigated. Let us distinguish the two following cases: 

(a) C+. In this case Vw = w(0+) G]0, 1[, there exists 7 C (w) e]0, 5/3] such that A 2 | 
if 7 I 7c (w). If 7 = 0, then A 2 > 0, V7 e]0,2]. 
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D + T\ 




(i) The only stable point is A+. (ii) The (unstable) equilibrium £ + is connected 

with C + and T> + with two separatrices spiraling 
through each other. 




i i j) 

(iii) No stable point exists: the flow is generically (iv) The only stable equilibrium point is B-, but 
forced to move away from E+. there are also generical solutions driving through 

A- and moving away from £_. 




FIGURE 3. Phase portrait sketches for the different situations summarized in 
Figure ||. 
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(b) C_. In this case Vw = u(0 ) g] — 1,-2/3 [, there exists 7c ( M ) G]0, 2] such that 
A 2 § Oif7 | 7c(«). Ifu(O-) G [-2/3,0], then A 2 > 0. 
In both cases, the expression of je(u(0)) is given by 7e(«) = ^fSf'- 

(4) T>±(-^ L ,+z J ,0±) and — z 7 , 0±) are four equilibria, where z 7 is chosen such 

that the points lie on the set \y 2 + z 2 = 3. These equilibria exist only if 7 < 5/3. 
When 7 = 5/3 we recover the case A± discussed above. If 7 < 5/3, the eigenvalues 
are given by A x = 0, A 2 = fgj > and A 3 = 2 J^±^^& , The 

corresponding eigenvectors are vi(0, 0, 1), v 2 is the tangent vector to the boundary 
of the set E± and v 3 is transversal to the same boundary. If u = u(0) G] — 1, — 2/3[, 
there exists a 7© (u) G]0, 2/3] such that A 3 = when 7 = 723 (w)- Otherwise, A 3 < 0. 
The expression for 70 (u) is 7© (it) = . 

(5) £±(- 

"a+v^1(o±)' ^) an( * ( _ Q +v ^I(o±)' _z o, 0) are other four equilibria, lying on the 
interior of £±. Let us distinguish the cases + and 0~: 

(a) On S + , the points are admissible if 7 < 5/3 (the limit case 7 = 5/3 the two 
points coincides with A+) in particular, Vu = u(0 + ) g]0, 1[, the equilibria on S + 
are admissible when 7 < 7c (w) (in the limit case jc(u) the two points coincide 
exactly with C + ). Note that y < 0. 

(b) On £_, Vw = u(0~) g] — 1, — 2/3[, the points are admissible if 7d(V) < 7 < 
7c («), the limit cases coinciding with P_ and C_. Notice that y > 0. 

(6) Finally, only in the case 7 = 4/3, there exists a set of equilibria given by the points 
(0, ±a/3, s), Vs. Note that, at level s = 0, that gives exactly the points V±. The eigen- 
values are given by X 1 = —4, A 2 = 1, A 3 = 0, with eigenvectors v 1; corresponding 
to a transversal direction to the boundary of the level set s = const., v 2 the tangent 
direction to the same boundary, and v 3 (0, 0, 1). 

2.2. Singularity formation. The analysis of the flow near the set 5fi allows to complete the 
study of late time behavior of the solution. 

Theorem 2.2. Let V £ £0 an d xq < 0. For almost every choice of initial data for the system 
0l.5aD - fll.5dD with k = and H(0) G [x^ , 0[, 3t s > such that the solution admits a right 



maximal extension to the set [0, t s [. The scalar field <p diverges in modulus ast — > t s , together 
with its velocity cj). 

Proof. First, we cast the system into the compact framework of problem V (see Definition 
|231 ). Using the analysis performed in last Section, we know that there are three possible cases 
in S + and five possible cases in XL, for the flow of the system at x = 0. By inspection 
of all possible cases, it is a lengthy but straightforward task to conclude that there cannot 
exist cycles contained in the set 5f2 and this means that each solution of problem V can only 
approach the stable equilibria of the system, which are the points A + and £L lying on the 
"unphysical" boundary Sf2. Let us observe also that A + is admissible only when 7 < 5/3. 

Translated into the original problem formalism with unknown functions depending on co- 
moving time t, the aforesaid means that there are two possible generical endstates for collaps- 
ing solutions: 

(1) either <p -> +00, <p -> +00, | -> -VQ and & -> 0, 
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(2) or (f> -»• -oo, -> -oo, 4 -»• and ^ -> 0. 
In both cases iJ diverges to — oo and it is easy to show that this happens as the comoving 



time approaches a finite value t s > 0, resulting in singularity formation. Indeed, from ( |1 .5b| ) 
(recall k = 0) and (pT6|), we have 

(2.11) H=-^(y 2 + 1Z 2 ). 

Since the quantity in round brackets approaches the positive value 6, there exists some to > 
such that H(t) < —2H(t) 2 and therefore, by standard comparison argument in ODE theory, 
the function H(t) is not extendable beyond some finite value t s . □ 

Remark 2.8. As seen before, the singularity forms because the scalar field energy \4>(t) 2 + 
V(<f>(t)) diverges. It is interesting to investigate the late time behavior of the fluid energy p(t). 



This can be easily inferred by inspection of (1.5d), recalling k = and 



P 1 / 3 7 , 

- + a 



p \y 

In case (Q), the quantity into round brackets approaches the value A/y/6 and therefore p(t) ~ 
e -44/Ve as ^ _ 00? so the energy density diverges. Instead, in case (jl]) (which can happen 
only if 7 < 5/3) p(t) ~ e <K&i-§/^ f which implies that the energy density of the fluid 
diverges only if 7 > 2/3 and is bounded otherwise. However, since tends to 0, the energy 
density of the fluid produces a lower order curvature divergence with respect to the energy of 
the scalar field, and the latter dominates in the late time of the collapse. 

3. GRAVITATIONAL COLLAPSE MODELS 

The model studied attracts a great deal of attention in cosmology, since it can be seen as the 
study of backwards-in-time singularity for a FRW spacetime in an expanding phase - in other 
words, reversing time direction we can study the initial singularity of expanding cosmologies. 
Nevertheless, we can also consider these singular solutions as matter interior, performing a 
suitable match with an exterior spacetime, thereby obtaining models of collapsing objects, 
whose endstate can be investigated. 

The natural choice for the exterior is the so-called generalized Vaidya solution, which is 
essentially the spacetime generated by a radiating object (see [^Dp for a detailed description), 

(3.1) ds c 2 xt = " (l - 2M y' dU 2 -2dYdU + Y 2 dfi 2 , 

where M is an arbitrary (positive) function. The matching is performed along a hypersurface 
S = {r = r b }. The Israel junction conditions at the matching hypersurface with a FRW 
interior have been shown in [ |i0j |TT| l and amounts to require that Misner-Sharp mass M is 
continuous and ^ = on the junction hypersurface. Note that the equation of motion for 
the scalar field remains smooth on the matching hypersurface. 

The endstate of the collapse of these homogeneous perfect fluid stars is analyzed in the 
following theorem (compare with flTT] , Theorem 4.1]). 

Theorem 3.1. A homogeneous perfect fluid star generically collapses to a black hole. 
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Proof. The equation of the apparent horizon for the fiat FRW metric is given by r 2 a(t) 2 = 
1. If a is bounded, the junction surface r = r b can be chosen sufficiently small such that 
(l — ^M) is bounded away from zero near the singularity. As a consequence, the exterior 
portion of the spacetime admits null radial geodesies that can be extended in the past up to the 
singularity which is therefore naked; otherwise, if d 2 is unbounded, the trapped region forms 



and the collapse is covered into a black hole QlOl] . 
Now, using ( |1.5a| )-( |T35D together with (^o]), one has 



(3.2) 



a(t) a(t) 



l--( y 2 + 1 z 2 ) 



a(t) a(t) 

and since, as a consequence from Theorem [2T2| , the quantity into square brackets converges 
generically to the value —2, we can integrate Q3.2D with respect to t and using the fact that 
a(t) — > we conclude that a(t) generically diverges to — oo and hence the singularity is 
hidden inside a black hole for almost each choice of the initial data. □ 

Remark 3.2. For the models under exam, the formation of naked singularities is forbidden for 
almost every choice of the initial data. Nevertheless, there may be special situations where 
the apparent horizon does not form and the singularity is hence naked - simply, the theorem 
excludes that this happens generically. 

For instance, with reference to the discussion of the equilibria given in Subsection we 
know that there exists non generical solutions tending to the equilibrium C + (Figure |3(ii)| ), 
such that y 2 + ^z 2 — > 6w(0 + ). Using ( |2.11| ) and fl3.2| ) respectively we have that, if < 
u(0 + ) < y/3/3, the singularity forms in a finite amount of comoving time and is naked 
because d remains bounded. Since z — > goes to zero, this case basically corresponds to 
the special case in the situation where only the scalar field appears, and not the perfect fluid 
- indeed, compare with JTT| , eqns.(5.3)-(5.5)]). There may be also special solutions future 
asymptotic to T> + , but in this case, anyway, the quantity y 2 + ^z 2 is seen to converge to the 
value . Then ( [2.1 1[ ) ensures singularity formation in a finite time, but the collapse ends 
into a black hole, as can be seen using again (|3T2|). 

In summary, we can conclude that the behavior here exhibited is - for the same class of 
potentials - quite similar to the scalar field without matter studied in JTT|]. As well known, 
Oppenheimer-Snider model of homogeneous dust represents a classical case of covered sin- 
gularity and this feature can be removed introducing higher order terms in the mass func- 
tions that destroy homogeneity. We recall that dust spacetimes can be regularly matched 
with Schwarzschild solution and, in this context, ( p.l| ) can be seen as a generalization of 
Schwarzschild. Therefore, it can be said the models studied in the present paper can be con- 
sidered a scalar field generalization of the classical Oppenheimer-Snider solution, and the 
addition of a barotropic perfect fluid to the matter content has not affected the general picture, 
that basically remains the same of the classical massless case. 
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